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An automated and robust method for adiabatic population transfer and the preparation of an arbitrary quantum 
superposition state in atomic system using the oscillating dark states (ODS) is presented. Quantum state of a 
three-level A configuration atomic system oscillates periodically between two ground levels, when two pairs 
of classical detuning laser fields driving the system into the ODS under evolving adiabatic conditions. The 
decoherence of the ODS evolution is greatly suppressed, and the oscillation is very stable, therefore adiabatic 
population transfer and the preparation of an arbitrary quantum superposition state of atomic system can be 
completed accurately and conveniently. 
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I. INTRODUCTION 

Dark state or coherent population trapping (CPT), which 
is null eigenvalue state of the interaction Hamiltonian, has 
given rise to growing interests for past few decades, since 
its first discovery in 1976 lM- If a three-level A atomic 
system is coherently trapped in dark states, there is no ab- 
sorption of incident laser fields even in the presence of reso- 
nant transitions, due to destructive quantum interference be- 
tween two transition pathways. Well-known dark state or 

CPT iSHBB&Slill] involves extensive applications in 
many fields such as electromagnetically induced transparency 
(EIT) B10I1 . lasing without inversio n I 1 11 1 1 2, 1 1 311 . the enhance- 
ment of the refractive index Il4 Il5ll . adiabatic population 
transfer |[l6[|T7tl. subrecoil laser cooling [18], and atom inter- 
ferometry 11911 . Coherent population transfer usin g st imulated 
Raman adiabatic passage(STIRAP) technique ll3, Tl7il20ll2lll . 
make it robust and efficient to prepare atoms and molecules 
in a well defined required state in many fields such as spec- 
troscopy, collision dynamics, and atomic physics. The Stokes 
and pump laser fields in a kind of counterintuitive order un- 
der two-photon resonance and adiabatic evolution carry out 
complete coherent population transfer, and the STIRAP tech- 
nique is insensitive to variations of laser pulse shape, intensity, 
and laser frequency. It is also of crucial importance that a re- 
quired coherent superposition state is prepared for atoms and 
molecules. Fractional STIRAP technique 1E2II23I1 . where the 
Stokes laser arrives before the pump laser and both terminate 
simultaneously maintaining a constant ratio of relative ampli- 
tudes, can create an arbitrary required coherent superposition 
state. The technique is also insensitive to pulse delay, intensity 
and frequency of laser fields. Another adiabatic transfer 112411 
using a tripod linkage by three laser fields in a four-level sys- 
tem can also to obtain good result. 

Here we propose an automated method for adiabatic pop- 
ulation transfer and the preparation of an arbitrary quantum 
superposition state in a three-level A configuration atomic sys- 



tem, using the oscillating dark states (ODS), where probability 
amplitudes of the system state oscillate periodically between 
two ground levels. The state to be prepared can be retrieved at 
a predictable time, and this method is insensitive to the initial 
state. 

The paper is organized as follows. The next section con- 
tains the discussion of the oscillating dark states (ODS) and 
evolving adiabatic conditions, under which the system follows 
the stable ODS evolution. In Sec. Ill, we introduce the master 
equation for a practical incoherent atomic system, and simu- 
late the ODS evolution, which then be used for adiabatic pop- 
ulation transfer and the preparation of an arbitrary quantum 
superposition state. We give some discussion and summarize 
our conclusions in Sec. IV. 



n. OSCILLATING DARK STATES 

The three-level A configuration atomic system is shown in 
Fig. Q] The ODS is implemented with two pairs of classical 
detuning laser fields coupling two lower ground levels |1) , 
|2) to a single upper level |3), respectively. We assume that 
the |1)-|2) transition is always dipole-forbidden. The time- 
dependent interaction Hamiltonian in the interaction picture 
that describes the atom-laser coupling within the dipole and 
rotating wave approximation (RWA) in the rotating frame 
reads (fi-l), 
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A=|(Aj + A 2 ), and A'=±[(A + A 2 ) - (A 3 + A 4 )]. A,=w 3 i- 
o>i(/=l,2), and A ; =« 3 2-w,0'=3,4) are the detunings of two 
pairs of laser angular frequencies w,0=l,2) and a>j(j=3,4) 
from the corresponding atomic transitions <y 3 i and co^- Real 
£2, and 0; 0=1, 2; 3, 4) are corresponding Rabi frequencies 
and phases of two pairs of laser fields, and we assume that 

Q.\=Q.2=£l\2, Q 3 =Q4=Q 3 4, <p\=<p2+lt=4>\2, an d 03=04=034, 

with the difference A0=0i2-0 3 4. The eigenstates of the in- 
teraction Hamiltonian (Q]), which can be generally described 
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FIG. 1: Two pairs of classical laser fields Cli , Q. 2 and fl 3 , fS 4 couple 
two lower ground levels |1) , |2) to a single upper level |3) , respec- 
tively, with their corresponding detunings A[ , A 2 and A 3 , A 4 , and 
their fixed phases: 0i=0 2 +7r, 3 =0 4 , in a three-level A configuration 
atomic system. 



in terms of the "mixing angles"^ and tp , are given as 
\ a+ ) = -/e /A0 sin6'sin^[l)+cos6lsin^|2)-cos 1; oe"'' 014 |3) , (2) 



compared with the separation of the corresponding eigenval- 
ues jum, 

\y6\ 2 <^\A -A ± \ 2 . (8) 
at 

Let Qi2=£l?4=Q (fixing Rabi frequencies when the system 
evolves) and Ai-A2=A3-A 4 =2( l >, thus we can obtain evolving 
adiabatic conditions according to ([8]) as 

\S\ 2 <k -|A ± VA 2 +4Q 2 | 2 , (9) 

which changes into \6\ <kQ, when A = 0, i.e. Ai=-A2 and 
A3=-A4. Evolving adiabatic conditions for the ODS ensure 
that the values of all parameters can be established in advance, 
because the evolution process is automated by laser-atom sys- 
tem itself, in contrast to ordinary dark states Q5|, Ij/ZD ^ where we 
manually with experimental instruments change relative val- 
ues of two Rabi frequencies to attain different dark states and 
to utilize STIRAP technique. 



|fl > = cos(9|l> - /sin<9e-'' A0 |2> , (3) 

|a_> = -ie iM> sin 6 cos + cos 6>cos</>|2> + sin (pe~^ 3i \3) , 

(4) 

where the mixing angles that dependent upon time t and Rabi 
frequencies are written as 

Q 12 sin(^f) 
tanfl = jV, (5) 

tan lip = — , (6) 

when we meet the condition A' = 0, i.e. (Ai + A2)=(A3 + A4). 
The eigenstate \oq) has no contribution from |3), and this 
ODS corresponds the null eigenvalue of interaction Hamilto- 
nian (Q~|l. If a system is in state \do), there is no possibility of 
excitation to |3) and subsequent spontaneous emission, thus 
the atom is effectively decoupled from laser fields. The eigen- 
values of the pair of states \a ± ) which contain a component 
of all three bare atomic states, are shifted up and down by an 
amount A ± , 

A ± = i[A ± VA 2 +4|P| 2 + 4|2| 2 ] . (7) 

The probability amplitudes of the system state can oscillate 
periodically between two ground levels |1) and |2) upon time 
f, when the system is driven into the ODS by laser fields with 
relatively fixed detunings and phases of four laser fields. 

Evolving adiabatic conditions, under which the system 
avoids diabatic coupling to bright states and is always in some 
dark states if it is initially in a certain one, is considered, since 
the ODS follows the evolution of a time-dependent Hamilto- 
nian (fTJ. The change rate of the mixing angle 9 must be small 



III. ADIABATIC POPULATION TRANSFER AND 
QUANTUM STATE PREPARATION 

A process for adiabatic population transfer and quantum 
state preparation is as follows: first, the initial instantaneous 
eigenstate of the whole system, which is determined by cho- 
sen interaction Hamiltonian, must be the atomic initial state, 
and here we assume that the atomic initial state is Then 
we concurrently and rapidly upload two pairs of laser fields 
as shown in Fig.Q] and the system periodically evolves subse- 
quently. Rapid unloading all laser fields simultaneously com- 
pletes a perfect population transfer from 1 1) to |2) after (1+f) 
(n-Q, 1,2, integral number) oscillation periods(r=27r/<5) if 
we want to make an adiabatic population transfer, and re- 
trieves a required quantum superposition state after time of 
(fo+«r)(?o is evolution time when the system firstly evolves 
from 1 1 ) to our required state, which can be easily predicted 
from ^ where a fixed A(f> should be calculated in advance) if 
we use the same experimental setup for the preparation of an 
arbitrary quantum superposition state. 

When we upload/unload two pairs of laser fields, up- 
loading/unloading adiabatic conditions are necessary besides 
evolving adiabatic conditions (0, in order to keep the least 
diabatic coupling between bright states and dark states of the 
system. We assume a concurrent increase/decrease of four 
laser fields in a short time t(t<zzT), i.e. £2i =£22=^3=^4=^ 
and j t Q\ = j,^- J,^~ Jt^- Thus, we can obtain the 
same result as evolving adiabatic conditions in (0 for upload- 
ing/unloading adiabatic conditions, according to ([8]). 

Time-dependent interaction Hamiltonian in interaction pic- 
ture that describes atom-laser coupling can be written as 

H' im = -I[O ie -<V Al W +Hz«-*»« i4 "cr 3 , 
+fi 3 e^V A3 'o-32 + iV^V^'c^ + H.c] , (10) 
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FIG. 2: Probability amplitudes evolution of initial state |1) finds 
expression in the time-dependent populations of three levels p\\, 
P22, P33, and coherence term between two ground levels |p2i| 2 with 
Aj = A 3 = 0.3, A, = A 4 = 0.2, y 2deph = 0.02, T 21 = 0.002, for 
different case: (a). SI = 2, (b). SI = 0.2, (c). Si = 0.08, respectively. 
Randomly chosen A(f> (adopted for different required superposition 
state, but having no effect on coherent population transfer process) 
gives the same periodical oscillation. All parameters are in the units 
of 7 3 i, and total evolution time t=4T is marked. 



where <x, ; =|/)01 (z. /=1,2, 3) are atomic projection operators. 
A master equation 10, H^l, using density matrix p, 

dp 
dt 
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is required, because the decoherence of a practical atomic sys- 
tem must be considered. T 3 i and F 3 2 represent rates of sponta- 
neous emission from level |3) to |1) and |2), respectively, and 
rates j2deph an d 7Meph describe energy-conserving dephasing 
processes. F2i is longitudinal relaxation from level |2) to |1), 
and it is much small compared with yzdeph- We can define co- 
herence decay rates as y3i=r" 31 +r 32 +y 3rfep /„ yii^n+jTdeph, 
and estimate F2i at 1/10 of jideph in the following analysis. 

The evolution of the ODS is simulated numerically as 
shown in Fig. [2] when the initial system is in [ 1) and a fixed 
A(f> should be calculated in advance (for a required arbitrary 
final state), and probability amplitudes evolution of the ODS 
with different Rabi frequencies Q is analyzed. The popula- 
tions of two ground lower levels, p\ \ and p 2 2, oscillate period- 
ically after we applying interaction Hamiltonian fllOt . and the 
population of upper level p 33 stays almost zero in total evo- 
lution stage, when evolving adiabatic conditions (O are well 



satisfied for the case in Fig.[2fa). The same periodical oscilla- 
tion of the coherence term |p 2 i| 2 between two ground levels as 
pn and P22, and |p 2 i| 2 ~pirp 2 2, imply a coherent evolution of 
quantum superposition states, i.e. the ODS. In case (b), Rabi 
frequency Q is decreased, and evolving adiabatic conditions 
are not well met, which results in a slight occupation of pop- 
ulation p 33 , but the oscillations of population p\\, pn and of 
coherence term IP21I 2 still dominate the whole evolution pro- 
cess. When we decrease Q further, laser fields are too weak to 
make upper level |3) occupied despite the violation of evolv- 
ing adiabatic conditions, and the coherence term of density 
matrix IP21I 2 keeps all the time in case (c), which implies 
there is no coherent evolution of the system state. 

From case (c) to (a), the value of Rabi frequencies CI of 
laser fields increased little by little, it can be seen that coher- 
ent reconstruction mechanism of laser fields becomes superior 
gradually to the decoherence effect of the system in case (c) 
and (b), and entirely dominates periodical oscillation of the 
system in case (a). A complete population transfer between 
two lower ground levels |1) and |2), and the preparation of an 
arbitrary quantum superposition state in an automatical man- 
ner, using a coherent evolution process of the ODS in case 
(a), is reliable. Here upper level |3) population is much lit- 
tle and the oscillation of probability amplitudes of two ground 
levels |1) and |2) is periodical, since the decoherence between 
two ground levels |1) and |2) can be greatly suppressed, ow- 
ing to coherent reconstruction contribution from laser fields. 
It can be seen that there are many optional retrieval time, with 
retrieval interval T/2 for an adiabatic population transfer pro- 
cess and with interval T for the preparation of an arbitrary 
quantum superposition state. 

With the contribution from applied laser fields, the state of 
the system evolves periodically and stably despite decay fac- 
tor j2\ ■ t expressed in coherence term of the system density 
matrix |p2i| 2 - When two pairs of laser fields are applied, we 
can also consider the long time behavior of state evolution 
of the system and analyze the fidelity of initial state 1 1 ) and 
retrieved states of the system after different evolution time t 
in Fig. [3] where only integral period evolution (t-nT) is ob- 
served. 721 is fixed at 0.022 in the units of y 3 i (reference to 
D\ line of ^Rubidium withy 3 i~36.10x 10 6 s _1 and estimated 
T2<fep/!=0.02y 3 i), and different total evolution times t is given. 
It can be seen that, when total evolution time t is up to long 
enough 1000r, where y 2 i ■ f>2000, high fidelity above 0.95 
is surprisingly available, which implies almost a coherence 
preserved and pure periodical evolution of the system. Fi- 
delity F decreases slightly in the first one T, during which a 
certain complicated amplitude decay of the ODS experiences, 
and stays unchangeably after it, which also can be found in the 
very beginning of evolution process in Fig. [2] High fidelity F 
mainly depends on relatively large value of Q. and always can 
be obtained by controlling laser fields. 



IV. DISCUSSION AND CONCLUSION 

The stability of frequency difference 26, n phase difference, 
and A<p of each pair of laser fields determine mainly the preci- 
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FIG. 3: Fidelity of initial state |1) and retrieved states of the system 
after different evolution time of the system (only integral period evo- 
lution observed, and up to a maximum 1000T here), with Q = 2, 
A] = A 3 = 0.3, A, = A 4 = 0.2, y 2] = 0.022. After the first one T 
(see inset), High fidelity above 0.95 is available and keep stable at all 
times, even if y 2 1 • f» 1 , where different initial states of system almost 
give the same result. 



sion of the ODS evolution. The population transfer efficiency 
and an ideal required superposition state mainly depend on 
the exact retrieval time that can be easily calculated according 
to (01. Two pairs of laser fields should be uploaded/unloaded 
spontaneously instead of counterintuitive and delayed upload- 
ing of one pair Stokes and pump laser fields. A concurrent 
uploading/unloading time t~ 10~ & s (a reasonable 0.01 T in our 
analysis which is short enough to avoid a certain decay of 
p at this stage) can be carried out using an optical switch, 
provided all the parameters are chose as ours, which can be 
tuned to change the oscillation period T. Laser fields may not 
be monochromatic, but the stability of relative phases of two 
pairs of laser fields, which can be easily achieved by acousto- 
optical modulation(AOM) device, can also guarantee almost 
the same result, despite a slow collective phase drift. A cho- 
sen 6 whose value is larger than the linewidth of laser fields 
can avoid undesirable interference of laser fields and optical 
pumping. The ODS method is immune to Doppler frequency 
shift, pulse shape of laser fields, and concurrent fluctuation of 
Rabi frequencies £2, 0=1, 2; 3,4). Well developed laser tech- 
niques can meet all above requirements to obtain a perfect adi- 
abatic population transfer and the preparation of an arbitrary 
quantum superposition state in atomic system. The greatly 
suppressed decoherence of the system according to our nu- 



merical calculation suggests the stability of the ODS evolu- 
tion. Numerous optional retrieval time owing to the property 
of the ODS evolution is available, with retrieval interval T/2 
for an adiabatic population transfer process and with interval 
T for the preparation of an arbitrary quantum superposition 
state, respectively. The transfer process is insensitive to the 
initial state of the atomic system 12711 . i.e., we can carry out 
complete coherent population transfer and obtain a required 
state at right retrieval time based on precise ODS evolution of 
atom-laser system, whatever the initial state of atomic system 
is, which is totally different from the STIRAP and fractional 
STIRAP technique. Different initial states of the system be- 
sides |1) here, to which the corresponding initial Hamiltonian 
of the system can be chosen to match the initial instantaneous 
eigenstate, give almost the same good result. Optical pumping 
is not necessary any longer in contrast to conventional STI- 
RAP technique, since it is not essential that initial state of 
the system must be any single lower state |1) or |2) when the 
automated ODS method is adopted. The ODS can be imple- 
mented in a three-level atomic system by using at least three 
laser fields in an appropriate configuration, but we use two 
pairs of laser fields in a symmetric scheme, in order to ob- 
tain complete oscillation of probability amplitudes between 
two ground levels. The ODS is similar to the famous Rabi 
resonance in a two-level atomic system driven by one reso- 
nant laser field @]i where we can not find the same preserved 
coherence. 

In conclusion, we propose an automated method using the 
ODS for a perfect adiabatic population transfer and the prepa- 
ration of an arbitrary quantum superposition state in atomic 
system. The preserved coherence, the automated manner, nu- 
merous optional retrieval time, an arbitrary initial state of sys- 
tem, make the ODS a potential method for physics process in 
atomic physics. 
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